Biquadratic interactions and spin-Peierls transition in the spin 1 chain LiVGe206 
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The magnetic susceptibility of a new one-dimensional, S=l system, the vanadium oxide LiVGe206, 
has been measured. Contrary to previous S=l chains, it exhibits an abrupt drop at 22 K typical of 
a spin-Peierls transition, and it is consistent with a gapless spectrum above this temperature. We 
propose that this behaviour is due to the presence of a significant biquadratic exchange interaction, a 
suggestion supported by quantum chemistry calculations that take into account the quasi-degeneracy 
of the t2g levels. 



The physics of one-dimensional (ID), spin 1 chains has 
attracted a considerable amount of attention after the 
prediction by Haldane that the ID Heisenberg model has 
a spin gap for integer spins a prediction confirmed 
since then by the observation of a gap in many spin 1 
chains A lot of progress has been made by study- 
ing the most general Hamiltonian describing an isotropic 
coupling between neighbouring spins 1, namely 

H^J'Y,SJ,+l+J"Y,{S^■S,+ lf (1) 
i i 

The phase diagram of this model is extremely rich [Q. 
It is most easily described using the parametrization 
J' = J cos 6, J" = JsinO. The system has two gapped 
phases: The Haldane phase for — 7r/4 < 9 < tt/A and a 
dimerized, gapped phase for — 37r/4 < 9 < — 7r/4. The 
two phases are connected by a critical point at 9 = — 7r/4 
for which the model is solvable by Bethe ansatz and 
the spectrum gapless. In the following, we will use the 
other standard parametrization J' = J, J" — —(5 J be- 
cause /3 is a direct measure of the relative strength of 
the biquadratic coupling to the bilinear one. With this 
parametrization, the Haldane phase is limited to J > 
and — 1 < /3 < 1, while the dimerized, gapped phase cor- 
responds to /3 > 1 if J > and /3 < — 1 if J < 0. Hope- 
lessly, it has not been possible so far to explore this phase 
diagram experimentally due to the lack of systems with 
a sizable biquadratic exchange: A well admitted point 
of view is that /3 is always very small, and the quantita- 
tive interpretation of the experimental results obtained 
on many spin 1 chain compounds with the pure Heisen- 
berg model has largely confirmed this assumption. 

In this Letter, we study the magnetic properties of a 
compound which has forced us to abandon this point of 
view, namely the vanadium oxide LiVGe206. Although 



this system is clearly a spin 1 chain given its structure, 
its susceptibility is markedly different from that of other 
spin 1 chains, and, as we shall explain below, the most 
plausible explanation seems to be that there is a signifi- 
cant biquadratic exchange. 




FIG. 1. Polyhedra representation of LiVGe206 structure 
down the a* axis. 

A powder sample of LiVGe206 was prepared by solid 
state reaction method starting from a stoichiometric mix- 
ture of Li2C03, Ge02 and V2O3 which was placed in a 
platinum crucible and heated under vacuum 2 hours at 
200° C and then 12 hours at 900° C. A second ther- 
mal treatment at 800° C for 12 hours was necessary to 
obtain the pure compound. Lithium vanadium metager- 
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manate LiVGe206 is the end member of the sohd solu- 
tion LiV(Si2-3;Gea;)06 (0 < x < 2) isolated and char- 
acterized by Satto et al ||]. This compound belongs 
to the alkali metal pyroxene family denoted AM'^+B2 06 
(A=alkali, B=Si, Ge, M=a variety of cations at a valence 
state 3-I-) which has been extensively studied since many 
compounds are naturally occuring minerals. The struc- 
ture of LiVGe206 has been determined from X-ray pow- 
der diffraction Rietveld analysis at 295 K. It crystallises 
in the monoclinic system, space group P2i/c compared 
to C2/c for LiVSi206 with a=9.863(4)A, b=8.763(2)A, 
c=5.409(l)Aand /3=108.24(1)°. The structure, depicted 
in Fig. 1 , is made of infinite isolated chains of edge shar- 
ing VOe octahedra linked together by chains of corner 
sharing Ge04 tetrahedra. The consequence of the change 
in the space group is that the chains of tetrahedra equiv- 
alent in LiVSi206 are distorted in LiVGe206 with two 
distinct crystallographic sites Gel and Ge2. The geo- 
metric parameters are reported in Table 1. Note that 
the chains are only connected to their neighbours through 
two Ge04 tetrahedra so that the coupling perpendicular 
to the chains must be very small. 
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FIG. 2. Magnetic susceptibility of LiVGe206. 

The magnetic susceptibility has been measured on a 
powder sample between 1.75 K and 350 K and for mag- 
netic fields ranging from 1 Oe to 7 T using a SQUID 
magnetometer. The results are insensitive to the mag- 
netic field strength, and the data are shown in Fig. 2. 
Although the impurity contribution at low temperature 
is not negligible, two features are clearly visible. The 
most spectacular one is an abrupt drop below 22 K. The 
behaviour below that temperature is typical of a spin- 
Peierls system with impurities [||, and it is very likely 
that this is indeed a spin-Peierls transition. In particu- 
lar, any explanation based on level crossings induced by 
the magnetic field can be excluded since the results are 
independent of the strength of the magnetic field. Be- 



sides, the susceptibility above this temperature is typical 
of gapless systems like the Heisenberg S=l/2 chain, and 
not of spin 1 chains with a Haldane gap, in which case 
the susceptibility drops much faster below its maximum. 
In fact, since the spin-Peierls instability occurs in gap- 
less systems because of the magnetic energy to be gained 
by opening a gap, both features are consistent with each 
other, and the main question is to understand the origin 
of the absence of a Haldane gap in the first place. 

The first possibility is that although vanadium is in 
the oxydation state 3-1- , i.e. the electronic configura- 
tion is 3cf, one could have an effective spin 1/2 model, 
as was advocated a long time ago by Castellani et al 
[0 in the context of V2O3. To check this, we have 
fitted the high temperature data with the expression 
x{T) — xvv + C /{T + 9). A reasonable fit is obtained 
for xvv = 6.3 X 10"'' cmVmol, C = 0, 80 cm^K/mol 
and 9 = 79K. This value of C is consistent with a spin 1 
system with an average g value of 1.79, a value typical of 
spin 1 V'^^. Besides, we have performed an ESR investi- 
gation of this system. The only signal that we could de- 
tect is typical of , but this signal corresponds to the 
small, paramagnetic contributions at low temperature. 
The fact that the main signal has not been detected is 
consistent with a spin 1 state originating from V^^ since 
any single-ion anisotropy pushes the resonance frequency 
outside the experimental range unless one works in a suf- 
ficiently high magnetic field. So it is clear that this is 
indeed a spin 1 system. 
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FIG. 3. Schematic representation of two-neighbouring V"^"'' 
ions with their surrounding oxygen octahedra. The solid 
(dashed) lines correspond to short (long) bonds. 

Taking for granted that the main contribution comes 
from spin 1 originating from V'^^, the possible expla- 
nations for the absence of a gap are: i) The coupling 
between the chains; ii) A strong uniaxial, easy plane 
anisotropy; iii) A significant biquadratic interaction. To 
review these various possiblities, it is useful to start with 
a brief description of the quantum chemistry of the vana- 
dium ions. In this compound, the vanadium atoms are 
sitting in a distorted oxygen octahedron (see Fig. 3) with 
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two neighbouring VO bonds shorter than the other four. 
This situation is very similar to the distortion of the sul- 
fur octahedron surrounding vanadium in AgVP2S6 |Q, 
and the three t2g orbitals are split into a low-lying dou- 
blet {dxy, dyz) and a single orbital {dxz) at an energy A 
above the doublet. Ab-initio calculations of the energy 
levels have been performed for a single ion embedded in 
a crystal, and the estimates for A are in the range 270- 
300 meV depending on the embedding. The degeneracy 
between the dxz and dyz orbitals is actually lifted due 
to the low symmetry, but this splitting is very small ac- 
cording to ab-initio calculations (in the range 3-7 meV) 
and physically unimportant, so we will neglect it for sim- 
plicity. The other on-site parameters are the Hund's rule 
coupling Jh and the repulsion between electrons on the 
same (resp. different) orbitals C/i (resp. U2)- Typical val- 
ues are Jh = 0.5-1 eV, Ui = 4-6 eV and U2 = Ui-2Jh. 
Vanadium first neighours sit in edge-sharing octahedra, 
the short bonds being opposite to each other and not in- 
volved in the common edge (see Fig. 3). If we denote by 
tij the hopping integrals between Wannier states corre- 
sponding to the t2g orbitals with i, j = 1, 2 and 3 for d^y, 
dyz and dxz respectively, then ti2 = ^13 = by symme- 
try while ^11, t22 and ^33 involve direct overlaps between 
3d orbitals and can be estimated from standard tables 
§: ill = 0.3 eV and ^22 = ^33 = -0.1 eV. Finally t23 
will have a contribution from indirect hopping through 
the oxygen 2p orbitals and is expected to be the largest 
hopping parameter in the problem, presumably in the 
range ^23 = 0.5 — 1 cV. Let us now discuss the possible 
explanations for the absence of a gap. 

i) Interchain coupling: It is well known that a rela- 
tively modest interchain coupling can close the gap . 
However, to have a significant effect on the susceptibility 
in the range 22-55K, the interchain coupling should be in 
that range. Such a large ratio of the interchain to intra- 
chain coupling seems very unlikely given the very one- 
dimensional structure of the compound. Besides if the 
interchain coupling was of the order of 20 K or larger, one 
should find a transition to long range order in the same 
temperature range, not a spin-Peierls transition. So we 
believe that the interchain coupling plays no significant 
role in the present system. 

ii) Uniaxial anisotropy: Several theoretical studies 
have shown that the presence of an easy plane, uniax- 
ial anisotropy of the form Dj^ii'^i)^ "^^^ close the Hal- 
dane gap, the critical value being D = J li the coupling 
between the spins is purely Heisenberg and isotropic in 
spin space ||ll[]. Can D be of the order of 50 K in the 
present system? The answer is clearly negative on both 
experimental and theoretical grounds. Given the elec- 
tronic structure of V'^"'" in the present compound, the 
spin-orbit coupling \L.S gives rise to a uniaxial, easy- 
plane anisotropy of order A^/A and a reduction of the 
perpendicular g factor: g± = 2(1 — A/A). In the isolated 
ion V'^"'", A = 104 cm~^ , which can be considered as an 



upper bound to its value in a crystal due to covalency 
effects. Now the average g factor determined from the 
high temperature susceptibility is equal to 1.79. Assum- 
ing that g\\ has the typical value 1.9, this implies that 
.g_L = 1.73, i.e. A/A = .13. With A = 104 cm'^, this 
leads to 13 = 20 K, a value clearly too small to close the 
gap. Besides, the actual value of A, hence the value of 
Z3, is certainly smaller than that. A careful analysis with 
single crystal NMR would be very useful in getting more 
precise estimates. Theoretically, the large value of A also 
points to a very small value of Z), of the order of a few K. 
This value is typical of V'^+ in this kind of environment 
|]l2| . So the explanation in terms of uniaxial anisotropy 
can be discarded. 

iii) Biquadratic interaction: Thus, the only re- 
maining possibility to explain the absence of a Haldane 
gap is the presence of a significant biquadratic interac- 
tion. Since the presence of a biquadratic interaction of 
significant magnitude has never been reported in any sys- 
tem, the most important task is to see whether such an 
interaction could indeed occur in the present case. So 
we have considered the general problem of the effective 
coupling between two spins 1 in almost degenerate t2g 
orbitals. This turns out to be a rather complicated prob- 
lem with very different results in different situations. So 
we will only present a summary of the results relevant for 
LiVGc206. The details will appear in a forthcoming pa- 
per Since the Hilbert space has a total of 495 states, 
this problem can be easily diagonalized numerically. As 
long as the energy A is not too small, the low energy sec- 
tor contains 9 states organized in 3 levels with degeneracy 
1, 3 and 5, typical of an isotropic coupling between spins 
1. Now the splitting between these levels can be used 
to determine both the bilinear and biquadratic exchange 
integrals since the energy levels of two spins 1 coupled 
by the Hamiltonian J[Si.S2 — /3(S'i.S'2)^] are given by 
Eo = -J(2-h4/3), = -J(l-h/3) and E2 = J(l-/3) for 
Stot = 0,1 and 2 respectively. It turns out that in a large 
parameter range, and in particular for physically relevant 
parameters, /? can be larger than 1, i.e. the biquadratic 
interaction can be significant and even dominant. For 
instance, if we take A — 300 meV, l}\ — h eV, — .5 
eV, ill = 0.3 eV and t22 = ^33 = —0.1 eV, /3 is already 
equal to .7 if i23 = -5 eV and increases with t23 until it di- 
verges when the bilinear exchange vanishes and becomes 
ferromagnetic. 

This remarkable behaviour is best understood in the 
context of a perturbation expansion in the hopping inte- 
grals. In the non-perturbed ground state each site carries 
two electons in a triplet state built out of the two low- 
lying orbitals. Then the degeneracy is in general lifted at 
second order, as for kinetic exchange, and the coupling 
is of the Heisenberg form, i.e. purely bilinear. What is 
important however is that it contains both ferromagnetic 
and antiferromagnetic contributions: t\\ and t22 yield an 
AF contribution Jaf = (l/2)(i?i + ti2)(l/£^2 + l/£^3) > 
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0, with E2 = Ui + and E-i = Ui + 2Jh + A due 
to Pauli exclusion principle, while ^23 yields a ferromag- 
netic contibution Jp = (2/3)^^3 (l/£;4 - i-/Ei) < with 
El =1/2 + A and E4 = U2 + 2Jh + A due to the Hund's 
rule coupling that favours intermediate states with high 
spin. Then if these couplings are of the same order, one 
has to go to fourth order to calculate the exchange inte- 
grals. It turns out that at this order the coupling is no 
longer purely bilinear but has a biquadratic contribution 
with a negative sign, and the ratio of this contribution to 
the bilinear one diverges at the transition between anti- 
ferromagnetic and ferromagnetic bilinear exchanges. The 
explicit expressions are very long and will be given in a 
forthcoming paper. In general, (3 becomes substantially 
large when the coupling from the second order contri- 
bution, \Jaf + Jf\, is comparable to ^^/{UfA), with t 
the largest of the hopping integrals. The parameters we 
estimate for LiVGe2 06 fall in this range. So, whereas 
all other explanations could be discarded on the basis of 
elementary but robust arguments, it seems that the pres- 
ence of a significant biquadratic interaction is a natural 
consequence of the relatively small splitting of the t2g 
levels, i.e. of the quasi orbital degeneracy of the system. 

Note that the order of magnitude of the exchange in- 
tegral supports by itself this explanation: Without any 
competition between Jp and Jaf the exchange integral 
would be much larger, hke e.g. in AgVP2S6, where it is 
of the order of 400 K |Q . As a consequence of this rela- 
tively small value of the nearest neighbour coupling, the 
relative magnitude of the next-nearest neighbour inter- 
action, which has been implicitly neglected so far, might 
actually not be as small as what a naive guess based on 
the geometry suggests. While such a term will, if any- 
thing, favour a dimerized groundstate, it does not close 
the Haldane gap and cannot by itself account for the 
properties of the present system. 

The presence of a biquadratic interaction has several 
interesting consequences. First of all, it will modify the 
temperature dependence of the susceptibility, and the 
data above the spin-Peierls transition are indeed consis- 
tent with a large value of f3. In particular, the susceptibil- 
ity is quite flat below the maximum, which is consistent 
with a very small gap, i.e. with a value of /3 around 1 
Besides, the gapped phase that occurs beyond 13—1 
corresponds to a two-fold, dimerized ground state, very 
much like for the spin 1/2 chain with large enough sec- 
ond neighbour exchange, and a Peierls transition with a 
dimerization of the lattice is of course a very likely in- 
stability of the system. So an explanation in terms of a 
biquadratic interaction of the same order but larger than 
the bilinear one is at the present stage the simplest and 
most natural explanation of the very peculiar properties 
of LiVGe206. Further experimental investigations of this 
material to confirm the Peierls nature of the transition 
and to study its magnetic properties with and without 
magnetic field is likely to be an active field of research in 



the next few years. 
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TABLE I. Selected geometrical parameters (A) 



4 



